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Abstract. We review several (and provide new) results on the theory of mo- 
ments, sums of squares and basic semi-algebraic sets when convexity is present. 
In particular, we show that under convexity, the hierarchy of semidefinite re- 
laxations for polynomial optimization simplifies and has finite convergence, a 
highly desirable feature as convex problems are in principle easier to solve. In 
addition, if a basic semi-algebraic set K is convex but its defining polynomials 
are not, we provide two algebraic certificate of convexity which can be checked 
numerically. The second is simpler and holds if a sufficient (and almost nec- 
essary) condition is satisfied, it also provides a new condition for K to have 
semidefinite representation. For this we use (and extend) some of recent re- 
sults from the author and Helton and Nie \E\. Finally, we show that when 
restricting to a certain class of convex polynomials, the celebrated Jensen's 
inequality in convex analysis can be extended to linear functionals that are 
not necessarily probability measures. 



1. Introduction 

Motivation. This paper is a contribution to the new emerging field of convex 
semi-algebraic geometry, and its purpose is threefold: First we show that the mo- 
ment approach for global polynomial optimization proposed in [13] , and based on 
semidefinite programming (SDP), is consistent as it simplifies and/or has better 
convergence properties when solving convex problems. In other words, the SDP 
moment approach somehow "recognizes" convexity, a highly desirable feature for a 
general purpose method because, in principle, convex problems should be easier to 
solve. 

We next review some recent results (and provide a new one) on the represen- 
tation of convex basic semi-algebraic sets by linear matrix inequalities which show 
how convexity permits to derive relatively simple and explicit semidefinite repre- 
sentations. In doing so we also provide a certificate of convexity for K when its 
defining polynomials are not convex. 

Finally, we consider the important Jensen's inequality in convex analysis. When 
restricting its application to a class of convex polynomials, we provide an extension 
to a class of linear functionals that are not necessarily probability measures. 

To do so, we use (and sometimes extend) some recent results of the author [161117] 
and Helton and Nie [B] . We hope to convince the reader that convex semi-algebraic 
geometry is indeed a very specific subarea of real algebraic geometry which should 
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deserve more attention from both the optimization and real algebraic geometry 
research communities. 

Background. I. Relatively recent results in the theory of moments and its dual 
theory of positive polynomials have been proved useful in polynomial optimiza- 
tion as they provide the basis of a specific convergent numerical approximation 
scheme. Namely, one can define a hierarchy of semidefinite relaxations (in short 
SDP-relaxations) of the original optimization problem whose associated monotone 
sequence of optimal values converges to the global optimum. For a more detail 
account of this approach, the interested reader is referred to e.g. Lasserre [T51 fT^] . 
Parrilo [2T], Schweighofer [2S], and the many references therein. 

Remarkably, practice seems to reveal that convergence is often fast and even 
finite. However, the size of the SDP-relaxations grows rapidly with the rank in the 
hierarchy; typically the r-th SDP-relaxation in the hierarchy has 0(n 2r ) variables 
and semidefinite matrices of 0(n r ) sizes (where n is the number of variables in the 
original problem). On the other hand, it is well-known that a large class of convex 
optimization problems can be solved efficiently; see e.g. Ben Tal and Nemirovski f J. 
Therefore, as the SDP-based moment approach is dedicated to solving difficult non 
convex (most of the time NP-hard) problems, it should have the highly desirable 
feature to somehow recognize "easy" problems like convex ones. That is, when 
applied to such easy problems it should show some significant improvement or a 
particular nice behavior not necessarily valid in the general case. Notice that this 
is not the case of the LP-based moment- approach described in pT] [15] for which 
only asymptotic (and not finite) convergence occurs in general (and especially for 
convex problems), a rather annoying feature. However, for SDP-relaxations, some 
results of [T7] already show that indeed convexity helps as one provides specialized 
representation results for convex polynomials that are nonnegative on a basic semi- 
algebraic set. 

II. Next, in view of the potential of semidefinite programming techniques, an 
important issue is the characterization of convex sets that are semidefinite repre- 
scntable (in short called SDr sets). A SDr set K C K™ is the projection of a set 
defined by linear matrix inequalities (LMIs). That is, 

n s 

K := {x G W : 3 y G R s s.t. A Q + V a A t + V y 2 Bj y 0} 

for some real symmetric matrices (Ai, Bj) (and where A y stands for A is positive 
semidefinite). For more details, the interested reader is referred to Ben Tal and 
Nemirovski pQ, Lewis et al. [TS], Parrilo [22], and more recently, Chua and Tuncel 
[2J , Helton and Nie [6] [7] , Henrion [8] and Lasserre [16] . For compact basic semi- 
algebraic sets 

(1.1) K := {x G R" : 9j (x) > 0, j = l,...,m}, 

recent results of Helton and Nie [6] [7] and the author [16] provide sufficient condi- 
tions on the defining polynomials (gj) c R[X] for the convex hull co (K) (= K if K 
is convex) to be SDr. Again, an interesting issue is to analyze whether convexity 
of K (with or without concavity of the defining polynomials (gj)) provides some 
additional insights and/or simplifications. Another interesting issue is how to de- 
tect whether a basic semi-algebraic set K is convex, or equivalently, how to obtain 
an algebraic certificate of convexity of K from its defining polynomials (jjj ). By 
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certificate we mean a mathematical statement that obviously implies convexity of 
K, can be checked numerically and does not require infinitely many tests. So far, 
and to the best of our knowledge, such a certificate does not exist. 

III. The celebrated Jensen's inequality is an important result in convex analysis 
which states that E fJi (f(x)) > f(E^(x)) for a convex function / : 1" — > K and a 
probability measure /i with E^(x) < oo. A third goal of this paper is to analyze 
whether when restricted to a certain class of convex polynomials, Jensen's inequality 
can be extended to a class of linear functionals larger than the class of probability 
measures. 

Contribution. Concerning issue I: We first recall two previous results proved in 
[17] : (a) the cone of convex SOS is dense (for the ^i-norm of coefficients) in the cone 
of nonnegative convex polynomials, and (b) a convex Positivstellensatz for convex 
polynomials nonnegative on K (a specialization of Putinar's Positivstellensatz). We 
then analyze the role of convexity for the polynomial optimization problem 

(1.2) P: /* = min { f(x) : x G K } 

X 

with K as in (11. 1| . and show that indeed convexity helps and makes the SDP- 
relaxations more efficient. In particular, when K is convex and Slater's condition^ 
holds, by using some recent results of Helton and Nie [6], we show that 

(i) If the polynomials /, (— gj) are all convex and V 2 / is positive definite (and 
so / is strictly convex) on K, then the hierarchy of SDP-relaxations has finite 
convergence. 

(ii) If / and (— gj) are all SOS-convex (i.e. their Hessian is a SOS matrix poly- 
nomial), then P reduces to solving a single SDP whose index in the hierarchy is 
readily available. 

Concerning II: Under certain sufficient conditions on the (gj) (typically some 
second order positive curvature conditions) Helton and Nie [H [7J have proved that 
co (K) (or K if convex) has a semidefinite representation that uses Schmiidgen or 
Putinar SOS representation of polynomials positive on K; see [6lll7j. Yet, in general 
its dimension depends on an unknown degree parameter in Schmiidgen (or Putinar) 
SOS representation. Our contribution is to provide a new sufficient condition for 
existence of a SDr when K is compact with nonempty interior and its boundary 
satisfies some nondegeneracy assumption. It translates the geometric property of 
convexity of K into a SOS Putinar representation of some appropriate polynomial 
obtained from each gj. When satisfied, this representation provides an algebraic 
certificate of convexity for K and it is almost necessary in the sense that it always 
holds true when relaxed by an arbitrary e > 0. It also contains as special cases 
Helton and Nie [6] sufficient conditions of SOS-convexity or strict convexity on 9K 
of the — gj's, and leads to an explicit semidefinite representation of K. We also 
provide a more general algebraic certificate based on Stengle's Positivstellensatz, 
but more complex and heavy to implement and so not very practical. In practice 
both certificates are obtained by solving a semidefinite program. Therefore, because 
of unavoidable numerical inaccuracies, the certificate is valid only up to machine 
precision. 



^Slater's condition holds for K in HI. Il l if for some xq G K, gj(xo) > 0, j = 1, . . . , m. 
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Concerning III, we prove that when restricting its application to the subclass of 
SOS-convex polynomials, Jensen's inequality can be extended to all linear function- 
als L y (with L y (l) = 1) in the dual cone of SOS polynomials, hence not necessarily 
probability measures. 

Some of the results already obtained in 6, 16] and in the present paper strongly 
suggest that the class of SOS-convex polynomials introduced in Helton and Nie [Bj 
is particularly nice and should deserve more attention. 

2. Notation, definitions and preliminary results 

Let R[X] be the ring of real polynomials in the variables X — (Xi, . . . , X n ), and 
let £ 2 [A] C R[X] be the subset of sums of squares (SOS) polynomials. Denote 
K[X]d C R[X) be the set of polynomials of degree at most d, which forms a vector 
space of dimension s(d) = ("+ d ). If / G R[X] d , write f(X) = £ aeN „ f a X a in 
the usual canonical basis (X a ), and denote by f = (f a ) £ R s ( d ) its vector of 
coefficients. Also write ||/||i (= ||f||i := ^3 \fa\) the ^i-norm of/. Finally, denote 
by £ 2 LY](i C S 2 [X] the subset of SOS polynomials of degree at most 2d. 

We use the notation X for the variable of a polynomial X i— » f(X) and x when 
a; is a point of K n , as for instance in {x G M n : f(x) > 0}. 

Moment matrix. With y = (y a ) being a sequence indexed in the canonical basis 
(X a ) of R[X], let L y : R[X] -> R be the linear functional 

/ (=J2f<*x a ) -> L y (f) = J2-f*y a , 

a a 

and let Md(y) be the symmetric matrix with rows and columns indexed in the 
canonical basis (X a ), and defined by: 

M d (y)(a,0) := L y (X a+fj ) = y a+f} , a,(3 6 N d 

with := {a G N" : |a| (= «i) < <*}• 

Localizing matrix. Similarly, with y = (y a ) and g G R[X] written 

X^g(X) = J2 a^x\ 

7 GN™ 

let M d (gy) be the symmetric matrix with rows and columns indexed in the canon- 
ical basis (X a ), and defined by: 

M d (gy)(a,P) := L y (g(X) X a +?) = ]T 57 y Q+;3+7 , 

7 

for every a,fl G NJJ. 

Putinar Positivstellensatz. Let Q(g) C R[X] be the quadratic module generated 
by the polynomials (gj) C R[X], that is, 

(2.1) Q(g) := j^o + ^ 5, : fo) C £ 2 [A] 1 . 

Assumption 2.1. K C R" is a compact basic semi- algebraic set defined as in 1 1.1]) 
and the quadratic polynomial X i— > M — \\ X\\ 2 belongs to Q(g). 
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Assumption 12. II is not very restrictive. For instance, it holds if every gj is affinc 
(i.e., K is a convex polytope) or if the level set {x : gj(x) > 0} is compact for some 
j G {1, . . . , to}. In addition, if M — \\x\\ > for all x G K, then it suffices to add 
the redundant quadratic constraint M 2 — \\x\\ 2 > to the definition (jTTTJ) of K and 
Assumption 12.11 will hold true. 



Theorem 2.2 (Putinar's Positivstellcnsatz [24 J. Let Assumption \2.1\ hold. If 

f G M.[X] is (strictly) positive on K 7 then f G Q(g)- That is: 



(2.2) 



3=1 



for some SOS polynomials (o~j )cS 2 [X] . 

2.1. A hierarchy of semidefinite relaxations (SDP-relaxations). Let P be 

the optimization problem (| 1 . 2[) with K as in (jl.ip and let Tj = [(deg <?j)/2] , j = 
1, . . . , to. With / G M[X] and 2r > max[deg /, maxj 2rj], consider the hierarchy of 
semidefinite relaxations (Q r ) defined by: 



(2.3) 



Qr 



inf L y (f) 
y 

s.t 



M r (y) h 
M r _ rj (g jy ) h0, 
2/o = 1 



j = 1,...,TO 



with optimal value denoted by inf Q r . One says that Q r is solvable if it has an 
optimal solution (in which case one writes inf Q r = min Q r ). The dual of Q r reads 



(2.4) 



sup A 
s.t 



f - \ = an + ^o-j gj 
a 3 eY?[X], j = 0,1,. 



, TO 

degcr , degtjj + degg-, < 2r, j = 1, . . . ,m 
with optimal value denoted by sup Q* (or max Q* if the sup is attained) . 



By weak duality sup Q* < inf Q r for every r G N and under Assumption 12.11 
inf Q r | /* as r — > oo. For a more detailed account see e.g. |13j . 



2.2. Convexity and SOS-convexity. We first briefly recall basic facts on a mul- 
tivariate convex function. If C C K™ is a nonempty convex set, a function / : C — * M 
is convex on C if and only if 

f(Xx + (1 - X)y) < Xf(x) + (1 - X)f(y), V A G (0, 1), x, y G C. 

Similarly, / is strictly convex on C if and only if the above inequality is strict for 
every x,y G C , x ^ y, and all A G (0, 1). 

If C C M" is an open convex set and / is twice differentiable on C, then / is 
convex on C if and only if its Hessian V 2 / is positive semidefinite on C (denoted 
V 2 / h on C). Finally, if V 2 / is positive definite on C (denoted V 2 / y on C) 
then / is strictly convex on C . 
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SOS-convexity. Helton and Nie [BJ have introduced the following interesting sub- 
class of convex polynomials, called SOS-convex polynomials. 

Definition 2.3 (Helton and Nie 6J). A polynomial f £ R[A]2d is said to be SOS- 
convex if V 2 / is SOS, that is, V 2 / = LL T for some real matrix polynomial L £ 
I[I]" XS (for some s £ N). 

As noted in [BJ, an important feature of SOS-convexity is that it can be can be 
checked numerically by solving a SDP. They have also proved the following impor- 
tant property: 

Lemma 2.4 (Helton and Nie [BJ Lemma 7]). // a symmetric matrix polynomial 
P eR[X] rxr is SOS then for any u £ R'\ the double integral 

X^ F(X,u) := / f P(u + s(X -u))dsdt 
Jo Jo 

is also a symmetric SOS matrix polynomial in H[X] rxr . 
And also: 

Lemma 2.5 (Helton and Nie [BJ Lemma 8]). For a polynomial f £ M[X] and every 
x,ueW n : 

f{x) = f(u) + Vf(u) T (x-u) 

+ (x-u) T / V 2 f(u + s(x - u))dsdt{x - u). 
Jo Jo 

" v ' 

F(x.u) 

And so if f is SOS-convex and f(u) — 0, V/(tt) = 0, then f is a SOS polynomial. 

2.3. An extension of Jensen's inequality. Recall that if fi is a probability 
measure on R" with E^ix) < oo, Jensen's inequality states that if / £ Li{jj) and 
/ is convex, then 

EM^)) > f( E M), 

a very useful property in many applications. 

We now provide an extension of Jensen's inequality when one restricts its ap- 
plication to the class of SOS-convex polynomials. Namely, we may consider the 
linear functionals L y : M.[X]-2d ^ R in the dual cone of X^LYJ^, that is, vectors 
y = [Va) such that Md(y) h and yo = L y (l) = 1; hence y is not necessarily the 
(truncated) moment sequence of some probability measure (i. Crucial in the proof 
is Lemma \2. 41 of Helton and Nie. 

Theorem 2.6. Let f £ R[A]2d be SOS-convex, and lety = {y a )aeN% d satisfy yo = 1 
and Mrf(y) y 0. Then: 

(2-5) L y (f(X)) > f(L y (X)), 

where L y (X) = (L y (Xx), . . . , L y {X n )). 
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Proof. Let z € 1™ be fixed, arbitrary, and consider the polynomial X ^ f(X) 
f(z). Then, 

(2.6) f(X)-f(z) = (Vf(z),X - z) + ((X - z),F(X)(X - z)), 

with F : R™ — > R[A]™ xn being the matrix polynomial 



A^ F(A) := / (\\ 
Jo Jo 



f(z + s(X -z))dsdt. 

As / is SOS-convex, by Lemma 12. 41 F is a SOS matrix polynomial and so the 
polynomial X h-> A(X) := ((X - z),F(X)(X - z) is SOS, i.e., A G £ 2 [A]. Then 
applying L y to the polynomial X i— > /(X) — /(z) and using ()2.6|1 yields (recall that 

yo = i) 

Ly(f(X))-f(z) = (Vf(z),L y (X)-z) + L y (A(X)) 

> (Vf(z),L y {X) - z) [because L y (A(X)) > 0]. 

As z G R™ was arbitrary, taking z := L y (X) (— (L y (Xi), . . . , L y (X n )) yields the 
desired result. □ 

As a consequence we also get: 

Corollary 2.7. Let f be a convex univariate polynomial, g 6 M.[X] (and so f o g £ 
R[A]J. Let d := [(deg / o ff)/2], and Zei y = (j/a)aeN" be such that yo = 1 and 
-Md(y) >r 0. Then: 

(2-7) iy[/(5W)] > /(i y [. 9 (A)]). 

Proof. Again let z € R n be fixed, arbitrary, and consider the univariate polynomial 
Y i ► /(F) - /(z) so that (J2ll) holds. That is, 

f{Y) - f{z) = f\z) (Y-z)+ F(Y)(Y - z) 2 , 

with F : R — > R[y] being the univariate polynomial 



Y^ F(Y) := f /*/" 
Jo Jo 



(z + s(Y - z))dsdt. 



As / is convex, /" > 0, and so the univariate polynomial Y \— » .F(Y)(Y — z) 2 is 
nonnegative, and being univariate, is SOS. Therefore, with Y := g(X), 

f(g(X)) - f(z) = f\z) (g(X) z) + F(g(X))(g(X) z) 2 , 

and so 

iy[/(.9(A))]-/(z) = f'(z)(L y [g(X)]-z)+L y [F(g(X))(g(X)-z) 2 ] 
> f(z)(L y [g(X)]-z) 

and taking z := L y [g(X)} yields the desired result. □ 



Hence the class of SOS-convex polynomials has the very interesting property 
to extend Jensen's inequality to some linear functionals that are not necessarily 
coming from a probability measure. 
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3. Semidefinite RELAXATIONS IN THE CONVEX CASE 

3.1. A convex Positivstellensatz. Let K be as in (jl.ip and define Q c (g) C M[X] 
to be the set: 

(3.1) Q c (g) := J a + A, g 3 ■ : A 6 M" 1 ; ct e £ 2 [X], ctq convex 1 C Q(g). 

The set Q c (g) is a specialization of Q(g) in (|2 . If) to the convex case, in that the 
weights asociated with the g/s are nonnegative scalars, i.e., SOS polynomials of 
degree 0, and the SOS polynomial <7o is convex. In particular, every / 6 Q c {g) 
is nonnegative on K. Let .Fk C M.[X] be the convex cone of convex polynomials 
nonnegative on K. 



Theorem 3.1 (Lasserre [H]). Let K be as in fl.l\) . Slater's condition hold and gj 
be concave for every j = 1, . . . , m. 

Then with Q c (g) as in 13.1]) . the set Q c (g) H^k is dense in Tk_ for the l\-norm 
|| • ||i. In particular, if K = R™ (so that Jg. —: T is now the set of nonnegative 
convex polynomials), then T, 2 [X] n T is dense in T . 

Theorem 13 . 1 1 states that if / is convex and nonnegative on K (including the case 
K = R n ) then one may approximate / by a sequence {f er } C Q c (g) H with 
11/ — fer\\i — > as e — > (and r — * oo). For instance, with r := [(deg/)/2j + 1, 

X^f cr (X) := f + e(9 ro (X)+9 r (X)), with 

r n y2k 

(3.2) X^0 r (X) := 1 + r ^ r - 

k=l i=l 

for some r e ; see Lasserre [T7J for details. Observe that Theorem 13 . 1 1 provides / with 
a certificate of nonnegativity on K. Indeed, let x £ K be fixed arbitrary. Then as 
fer S <3e(g) one has fer(x) > 0. Letting e J, yields < lim e ^ fer{ x ) = /O^)- And 
as a; € K was arbitray, / > on K. 

Theorem [3J] is a convex (weak) version of Theorem l2.2l fPutinar's Positivstellen- 
satz) where one replaces the quadratic module Q(g) with its subset Q c (g)- We call 
it a weak version of Theorem [5121 because it invokes a density result (i.e. f er £ Q c {g) 
whereas / might not be an element of Q c (<?)). Notice that / is allowed to be non- 
negative (instead of strictly positive) on K and K need not be compact; recall that 
extending Theorem 12.21 to non compact basic semi-algebraic sets K and to polyno- 
mials / nonnegative on K is hopeless in general; see Scheiderer |26j . 



Corollary 3.2. Let K be as in {OJ), / € R[X] with f* := rmn x {f{x) : x £ K} 
and let d :— max[|~(deg /)/2] , maxj |~(deg <?j)/2] ]. Consider the simplified SDP- 
relaxation 

inf L y (f) 

(3.3) Q : { * 

L y (g 3 ) > o, i = l,...,m 

Vo = 1 
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id its dual 



(3.4) Q* 



sup 7 

7,fo, A 

m 

s.t. /-7 = cr + ^A i5:) 

3=1 

a G S 2 [X] d ; Aj ■> 0, j = l,...,m 

(a) If f — f* S Qc(g) then the SDF '-relaxation Q and its di/aZ Q* are exact. 

(b) I//, — <7j 6 R[JT] are convex, j = 1, . ..,m, and ify is an optimal solution of 
Q which satisfies 

(3.5) rankAf d (y) = rankM d _i(y), 

i/ien Q is exact and x* := (L y (Xi)) G K is a (global) minimizer of f on K. 

Proof (a) If / — /* G Q c (<?), i.e., if /-/* = <r + £™ ^ A, 5j , with ao 6 Y?[X] d and 
A G R™, the triplet (f*,a , A) is a feasible solution of Q* with value /*. Therefore, 
as supQ* < inf Q < /*, the SDP-relaxation Q and its dual Q* are exact. In fact, 
(/*, (To, A) is an optimal solution of Q*. 

(b) If y satisfies the rank condition (|3.5p then by the flat extension theorem 
of Curto and Fialkow [4], y is the (truncated) moment sequence of an atomic 
probability measure /i on R™, say \i — Yfk=l ^fc^(fc) with s = rankMd(y), < 
Afc < 1, J2k = 1, and S x ^) being the Dirac measure at x(k) G R n , k = 1, . . . , s. 
Let x* := J2k ^kx{k) = (L y (X t )) G R". Then /* > L y (f) and by convexity of /, 
L y(f) = J2k X kf(x(k)) > /(Efe A fca;(fc)) = f(x*). Similarly, by convexity of -g h 
< L y {g j ) = J2k X k9 3 (x(k)) < gj(J2k X kx(k)) = gj(x*), j = l,...,m. Therefore, 
x* G K and as f(x*) < /*, x* is a global minimizer of / on K. □ 

Notice that K in Corollary 13.21 need not be compact. Also, Corollary 13.2T b) 
has practical value because in general one does not know whether / — /* G Q c (g) 
(despite that in the convex case, / — /* G J~k and Q c {g) H is dense in ^k)- 
However, one may still solve Q and check whether the rank condition (|3.5[) is 
satisfied. If in solving Q r , the rank condition (|3.5p is not satisfied, then other 
sufficient conditions can be exploited as we next see. 

3.2. The SOS-convex case. Part (a) of the following result is already contained 
in Lasserre [T71 Cor. 2.5]. 



Theorem 3.3. Let K be as in il.l)) and Slater's condition hold. Let f G R[X] be 
such that f* := iitf x {f(x) : x G K} = f(x*) for some x* G K. If f is SOS-convex 
and —gj is SOS- convex for every j = 1, . . . ,m, then: 

(a) /-/* eQ c (g). 

(b) The simplified SDP-relaxation Q in \3.3)) and its dual \3.$ are exact and 
solvable. Ify is an optimal solution of Q then x* := (i y (Xj)) G K is a global 
minimizer of f on K. 

Proof, (a) is proved in [T71 Cor. 2.5]. (b) That Q is exact follows from (a) and 
Corollary I3.2f a). Hence it is solvable (e.g. take y to be the moment sequence 
associated with the Dirac measure at a global minimizer x* G K). So let y be an 
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optimal solution of Q, hence with /* = L y (f). As — gj is SOS-convex for every j, 
then by Theorem l2~6l < L y (gj) < gj(x*) with x* := (£ y (X»)) and so x* G K. 
Similarly, as / is SOS-convex, we also have /* = L y {f) > f(x*) which proves that 
f(x*) = f* and x* is a global minimizer of / on K. Finally, as by (a) / — /* G Q c (g) 
then Q* is exact and solvable. □ 

(Again notice that K in Theorem 13.31 need not be compact.) So the class of 
SOS-convex polynomials is particularly interesting. Not only Jensen's inequality 
can be extended to some linear functionals that are not coming from a probability 
measure, but one may also solve SOS-convex optimization problems P in (|1.2p 
(i.e. with / and K defined with SOS-convex polynomials) by solving the single 
semidefinite program (|3.3|) . 

Notice that a self-concordan10 logarithmic barrier function exists for (|3.3|) whereas 
the logarithmic barrier function with barrier parameter /x: 

m 

(3.6) xh-K/> M (x) := nf(x) -^2ln(-gj(x)), 

3=1 

associated with P, is not self-concordant in general. Therefore, despite (|3.3p in- 
volves additional variables (a lifting), solving (|3.3p via an interior point method 
might be more efficient than solving P by using the logarithmic barrier function 
(|3.6[) with no lifting. In addition, all SOS-convex polynomials nonnegative on K 
and which attain their minimum on K, belong to Q c (g), a very specific version of 
Putinar Positivstellensatz (as / is only nonnegative and K need not be compact). 

3.3. The strictly convex case. If / or some of the — g^s is not SOS-convex but 
V 2 / >- (so that / is strictly convex) and — gj is convex for every j = 1, . . . , m, 
then inspired by a nice argument from Helton and Nie [6] for existence of a semi- 
definite representation of convex sets, one obtains the following result. 

Theorem 3.4. Let K be as in 1 1.1]) and let Assumption \2. 1\ and Slater's condition 
hold. Assume that /, — gj G M[X] are convex, j = 1, . . . ,m, with V 2 / >- on K. 

Then the hierarchy of SDF '-relaxations defined in 12. 3]) has finite convergence. 
That is, f* — supQ* = inf Q r for some index r. In addition, Q r and Q* are 
solvable so that f* — maxQ* = minQ r . 

Proof. Let x* G K be a global minimizer (i.e. /* = f(x*)). As Slater's condition 
holds, there exists a vector of Karush-Kuhn- Tucker (KKT) multipliers A G such 
that the (convex) Lagrangian Lf G M[X] defined by 

m 

(3.7) X i * L f (X) := f(X) - /* - £ X i 

j'=i 

has a global minimum at x* G K, i.e., VI/(z*) = 0. In addition, Xjgj(x*) = for 
every j = 1, . . . ,m and Lf(x*) = 0. Then, by Lemma 12.51 

L f (X) = ((X — x*),F(X, x*)(X - x*)) 



2 The self-concordance property introduced in 1201 is fundamental in the design and efficiency 
of interior point methods for convex programming. 
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with 

F(X,x*) := (J V 2 L f (x* +s{X -x*))dsd?J . 

Next, let I n be the n x n identity matrix. As V 2 / V on K, continuity of the 
(strictly positive) smallest eigenvalue of V 2 / and compactness of K yield that 
V / y SI n on K, for some 5 > 0. Next, as —gj is convex for every j, and in view 
of the definition of L f , V 2 L/ h V 2 / h 5I n on K. Hence for every ^R", 
£, T F(x,x*)£ > 6 f*£, T £,dsdt = §£ T £, and so F(x,x*) h §/„ for every x G K. 
Therefore, by the matrix polynomial version of Putinar Positivstellensatz, 

m 

F(X,x*) = F (X)+Y,Fj(X)9j(X), 

.7=1 

for some real SOS matrix polynomials X i— > Fj(X) = Lj(X)Lj(X) T (for some 
apppropriate Lj G M[X]™ XW ), j = 0, ...,m. See Helton and Nie [6], Kojima and 
Maramatsu [10], Hoi and Scherer [IT] . But then 

X i ► ((X-ir*),^^,^*)^ - a;*)) = ^(X) e S 2 [X], j = 0, . ..,m 

and so 

m 

/(X)-/* = L/(X) + ^A ift (X) 

m 

= a (X) + J2(^+^(X))9j(X). 

3=1 

Let 2s be the maximum degree of the SOS polynomials (oj). Then (/*, {<jj + A^}) 
is a feasible solution of the SDP-relaxation Q* in (|2.4|) with r := s + max^rj. 
Therefore, as supQ* < inf Q r < /*, the SDP-relaxations Q r and Q* are exact, 
finite convergence occurs and Q* is solvable. But this also implies that Q r is 
solvable (take y to be the moment sequence of the Dirac measure 5 X * at any global 
minimizer x* G K). □ 

When compared to Theorem 13.31 for the SOS-convex case, in the strictly con- 
vex case the simplified SDP-relaxation Q in (|3.3p is not guaranteed to be exact. 
However, finite convergence still occurs for the SDP-relaxations (Q r ) in (|2.3[) . 

Remark 3.5. It is worth emphasizing that in general, the hierarchy of LP-relaxations 
(as opposed to SDP-relaxations) defined in [TS] and based on Krivine's represen- 
tation [121 130] for polynomials positive on K, cannot have finite convergence, es- 
pecially in the convex case! For more details, the interested reader is referred to 
[14] IT5] . Therefore, and despite LP software packages can solve LP problems of 
very large size, using LP-relaxations does not seem a good idea even for solving a 
convex polynomial optimization problem. 

4. Convexity and semidefinite representation of convex sets 

We now consider the semidefinite representation of convex sets. First recall the 
following result. 
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Theorem 4.1 (Lasserre |16j). Let K in (HHj) be compact with gj concave, j = 
l,...,m, and assume that Slater's condition holds. If the Lagrangian polynomial 
Lf in |ff.7| ) associated with every linear polynomial f £ M.[X] is SOS, then with 
d := maXjf(deggj)/2], the set 

( M d {y) >Q 
< RS (2d) . I L y (gj) > 0, j = l,...,m 
| Ly(Xi) =a;,-, i = l, ...,n 
I y =1 

is a semidefinite representation o/K. 



(4.1) 



ft := ^ 



0,y) £ 



Next, Helton and Nie [6l [7] have provided several interesting second-order pos- 
itive curvature (sufficient and necessary) conditions on the defining polynomials 
(gj) for K (or its convex hull co(K)) to have a SDr. In particular (recall that 
r j = r( dc g5j )/ 2 l for every j = 1, . . . , m): 



Theorem 4.2 (Helton and Nie [6j). Let K in U.l\) be convex, Asssumption \2.l\ 
hold, and assume that Slater's condition holds and gj is concave onK, j = 1, . . . , to. 

(a) If —gj is SOS-convex for every j = 1, . . . , m, then for every linear f £ K[X], 
the associated Lagrangian Lf {3. 7| ) is SOS and the set Q in is a semidefinite 
representation o/K. 

(b) If every —gj is either SOS-convex or satisfies —V 2 gi >~ on Kfl{i : gj(x) — 
0}, then there exists r £ N such that the set 

" M r (y) h0 



(4.2) n.{M 6 rxR^) : I = ^ i = 1> n 

is a semidefinite representation o/K. 



M r- rj (g 3 y) to, j' = 1,...,to 
yo = f 



See [SJ Theor. 6, and 9]. This follows from the fact that the Hessian V 2 L/ 
associated with a linear / £ M[X] has a Putinar representation in terms of SOS 
matrix polynomials, and with degree of the weights bounded uniformly in /. In 
principle, the degree parameter d in Theorem 14. 2f b) may be computed by solving 
a hierarchy of semidefinite programs. Some other (more technical) weaker second- 
order positive curvature sufficient conditions (merely for existence of a SDr) are also 
provided in [6l [7] but the semidefinite representation is not explicit any more in 
terms of the defining polynomials (gj ) . Notice that if K is compact but Assumption 
12.11 does not hold, then one still obtains a semidefinite representation for K but 
more complicated as it is now based on Schmudgen's representation [27 instead of 
Putinar's representation; see [6l Theor. 5], 

We next provide a sufficient condition in the case where K is convex but its 
defining polynomials (—gj) are not necessarily convex. Among its distinguishing 
features, it is checkable numerically, contains Theorem 14.21 as a special case and 
leads to the explicit semidefinite representation (|4.2[) of K. 

4.1. Algebraic certificate of convexity. We first present the following charac- 
terization of convexity when K is closed, satisfies a nondegeneracy assumption on 
its boundary, and Slater's condition holds. 
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Lemma 4.3. Let K be as in fUZP (hence closed), Slater's condition hold and 
assume that for every j — 1, . . . , to, Vgj(y) ^ if y £ K and gj(y) = 0. Then K 
is convex if and only if for every j = 1, . . . ,m, 

(4.3) (V 9j (y),x-y) >0, Vx andVy eK with g 3 (y) = Q. 

Proof. The only if part is obvious. Indeed if (ygj(y),x — y) < for some a; G K 
and y £ K with <?j(y) = 0, then there is some t > such that gj(y + i(x — y)) < 
for all t £ (0, t) and so the point x' := tx + (1 — t)y does not belong to K, which in 
turn implies that K is not convex. 

For the if part, (|4.3|) implies that at every point of the boundary, there exists 
a supporting hyperplane for K. As K is closed with nonempty interior, the result 
follows from ^ Theor. 1.3.30. □ 

The nondegeneracy assumption is crucial as demonstrated in the following simple 
example kindly provided by an anonymous referee: 

Example 1. Consider the non convex set K C R 2 denned by: 

K := { x € R 2 : (1 - x\ + x 2 2 f > 0, 10 - x\ - x\ > } 

Then it is straightforward to see that (|4.3|l is satisfied. This is because Vgi vanishes on 
the piece of boundary determined by gi(x) — 0. 

Next, using the above characterization (|4.3[) . we provide an algebraic certificate 
of convexity. 

Corollary 4.4 (Algebraic certificate of convexity). Let K be as in il.l)) . Slater's 
condition hold and assume that for every j — 1, . . . , to, Vgj{y) / t/ j € K and 
9j(u) = 0- Then K is convex if and only if for every j — 1, . . . , m, 

(4.4) h 3 {X,Y){Vg 3 {Y),X-Y) = {V 9j (Y), X-Y) 21 +9j(X, Y)+<pj(X, Y) 9j (Y), 

for some integer I £ N, some polynomial tpj G WL[X, Y] and some polynomials hj, 9j 
in the preordering of R[X,Y] generated bythefamily of polynomials (gk(X), g p (Y)), 
k,p e {1, . . . ,m}, p ^ j. 

Proof. By Lemma 14 . 31 K is convex if and only if for every j = l,...,m, the 
polynomial (X,Y) ^ (Vgj(Y), X — Y) is nonnegative on the set VLj defined by: 

(4.5) % := {(x,y)eKxK : 9j (y)=0}. 

Equivalently, K is convex if and only if for every j — 1, . . . ,m: 

« = {Mel": (^)GKxK; 9j (y) = 0; 

{V 9j (y),x-y)<0 ] (V 9j (y),x-y)^0}. 

Then gl} follows from Stengle's Positivstellensatz [25l Theor. 4.4.2, p. 92]. □ 



^The author is grateful to L. Tuncel for providing us with the reference |28j . 
4 The preordering of R[X] generated by a family [g\, . . . , g m ) C R[X] is the set of polynomials 
{P ■■ P = £jc{i,.... m }*j(IW<&). witha /6 E 2 [V]}. 
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Observe that Corollary 14.41 provides an algebraic certificate of convexity when 
K is closed with nonempty interior and a nondegeneracy assumption holds on its 
boundary. If one fixes an a priory bound s on I g N and on the degree of hj , 9j and 
ifj, then checking whether (|4.4[) holds reduces to solving a semidefinite program. If 
K is convex, by increasing s, eventually one would obtain such a certificate if one 
could solve semidefinite programs exactly. In practice, and because of unavoidable 
numerical inaccuracies, one only obtains a numerical approximation of the optimal 
value and so, a certificate valid up to machine precision only. 

However, implementing such a procedure is extremely costly because one has 
potentially 2 x 2 m unknown SOS polynomials to define hj and 9j in (|4.4|) ! Therefore, 
it is highly desirable to provide a less costly certificate but with no guarantee to 
hold for every K as in Corollary 14.41 

In particular one only considers compact sets K. Indeed, if K is compact, one 
has the following result (recall that go = 1). 

Lemma 4.5. LetK. be convex, Assunwtion \2. 1\ and Slater's condition hold. Assume 
that for every j — 1, . . . , to, Vgj(y) /Oi/j/eK and gj(y) = 0. Then for every 
e > and every j = 1, . . . , to: 

m m 

(V gj (Y),X-Y)+e = ^a jk (X,Y)g k (X) + £ 1> Jk (X,Y)g k (Y) 

k=0 k=0,k^tj 

(4.6) +ip j (X,Y)g j (Y), 

for some SOS polynomials (ojk) and (ipjk)k^j C S 2 [A, Y], and some polynomial 
ipj €R[X,Y]. 

Proof. By Lemma |4.3[ for every j = 1, . . . , to, and every i,i/eK such that gj(y) = 
0, ()4.3p holds and therefore, for every j = 1, . . . , m, 

(4.7) (Vg 3 (y),x-y) + e > V(x,y) G %, 

where ilj has been defined in (|4.5| . As K satisfies Assumption 12 . II then so does fl,- 
for every j — 1, . . . , to. Hence (|4.6[) follows from (|4.7[) and Theorem 12.21 □ 

Therefore, inspired by Lemma 14.51 introduce the following condition: 

Assumption 4.6 (Certificate of convexity). For every j — l,...,m, j4-6\ ) holds 
with e = 0. Then let dj G N be such that 2dj is larger than the maximum degree of 
the polynomials Uj k g k , 4>jkgk, ?Pj9j <= R[-X"i Y] in \4-.0j ), j = 1, . . . , m. 

When K is closed (and not necessarily compact), Slater's condition holds and 
the nondegeneracy assumption on the boundary holds (i.e., ^7gj(y) ^ if y £ 
K and gj (y) = 0) Assumption 14.61 is indeed a certificate of convexity because 
then (|4.3[) holds for every x,y € K with gj(y) — 0, and by Lemma K is 
convex. It translates the geometric property of convexity of K into an algebraic SOS 
Putinar representation of the polynomial (X,Y) i— > (Vgj(Y),X — Y) nonnegative 
on Qj, j = 1 , . . . , to. On the other hand, if K is convex and Assumption 12. l\ 
Slater's condition and the nondegeneracy assumption all hold, then Assumption 
14.61 is almost necessary as, by Lemma l4~5| (|4.6j) holds with e > arbitrary. 
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With dj fixed a priori, checking whether (|4.6|) hold with e = can be done 
numerically. (However, again it provides a certificate of convexity valid up to ma- 
chine precision only.) For instance, for every j = 1, . . . ,m, it suffices to solve the 
semidefinite program (recall that r^ = |~(deg<7fc)/2], k — 1 . . . , m) 

min L x ((V gj (Y),X-Y)) 
M d] (z) t 

M d ^ rk (g k (X)z) bO, k = l,...,m 
M d ^ rk (g k (Y)z)hO, k = l,...,m;k^j 
M d ^ rj ( 9j (Y)z) = 
2/0 = 1 

If = for every j — 1, . . . , m, then Assumption 14.61 holds. This is in contrast 
to the PP-BDR property in [17] that cannot be checked numerically as it involves 
infinitely many linear polynomials /. 

Remark 4.7. Observe that the usual rank condition (|3.5p used as a stopping 
criterion to detect whether (|4.8p is exact (i.e. pi = 0), cannot be satisfied in 
solving (|4.8p with primal dual interior point methods (as in the SDP-solvers used by 
GloptiPoly) because one tries to find an optimal solution z* in the relative interior 
of the feasible set of (|4.8p and this gives maximum rank to the moment matrix 
M dj (z*). Therefore, in the context of (|4.8|1 . if indeed pj — then z* corresponds 
to the moment vector of some probability measure /i supported on the set of points 
(x, x) e K x K that satisfy gj(x) = (as indeed L z . ((Vgj(Y),X - Y))) = = pj). 
Therefore pj = as dj increases but the rank of M dj (z*) does not stabilize because 
p is not finitely supported. In particular, a good candidate z* for optimal solution 
is the moment vector of the probability measure uniformly distributed on the set 
{{x,x) E K x K : gj {x) = 0}. 

Alternatively, if pj w and the dual of (14. 8|) has an optimal solution {ojh ,^jki ipj ) , 
then in some cases one may check if (|4.6[) holds exactly after appropriate rounding 
of coefficients of the solution. But in general, obtaining an exact certificate (i.e., 
Pj = in the primal or (|4.6p with e = in the dual) numerically is hopeless. 

Example 2. Consider the following simple illustrative example in R 2 : 

(4.9) K := {x £ R 2 : Xl x 2 - 1/4 > 0; 0.5 - (xj - 0.5) 2 - {x 2 - 0.5) 2 > 0} 

Obviously K is convex but its defining polynomial x gi(x) := x\x 2 — 1/4 is not concave 
whereas x \— > g 2 (x) := 0.5 — (n — 0.5) 2 — (x 2 — 0.5) 2 is. 

With di = 3, solving f|4.8[) using GloptiPoly 3^ yield s the optimal value pi ~ —4.58.10 
which, in view of the machine precision for the SDP solvers used in GloptiPoly, could be 
considered to be zero, but of course with no guarantee. However, and according to Remark 
14.71 we could check that (again up to machine precision) for every a £ N" with \a\ < 2dj, 
z a,a = z 2a ,o and = z o,a- I n addition, because of symmetry, z a ,/3 = z a i whenever 
a[ — a 2 and a' 2 = ot\ (and similarly for (3 and /?'). Indeed for moments of order 1 we have 
<,3 = (0.5707, 0.5707, 0.5707, 0.5707) and for moments of order 2, 

z* a ,fi = (0.4090, 0.25, 0.4090, 0.25, 0.4090, 0.25, 0.4090, 0.4090, 0.25, 0.4090). 



(4.8) 



Pj ■ 
s.t. 



5 GloptiPoly 3 (a Matlab based public software) is an extension of GloptiPoly [9] 
to solve the generalized problem of moments described in |18| . For more details see 
www. laas . fr/~henr ion/software/. 
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For j = 2 there is no test to perform because — g 2 being quadratic and convex yields 
(4.10) (Vg 2 (Y),X-Y) = g 2 (X) - g 2 (Y) + (X - Y) T '(-V 2 g 2 (Y))(X - Y) 

SOS 

which is in the form (|4.6p with d 2 = 1. 

We next show the role of Assumption 14.61 in obtaining a semidefinite represen- 
tation of K. 



Theorem 4.8. Let Assumption \2. 1\ and Slater's condition hold. Moreover, assume 
that for every j = l,...,m, V<?j(y) ^ whenever y 6 K and (^(y) = 0. // 
Assumption \4-6\ holds then K is convex and £1 in with d := maxj dj, is a 

semidefinite representation o/K. 

Proof. That K is convex follows from Lemma 14.31 We next prove that the PP- 
BDR property defined in Lasserre [TB] holds for K. Let / G M.[X] be a linear 
polynomial with coefficient vector f £ K™ (i.e., X i— > /(X) = f T X) and consider 
the optimization problem P : min {{ T x : x G K}. As K is compact, let x* G K 
be a global minimizcr of /. The Fritz- John optimality conditions state that there 
exists ^ A G R" l+1 such that 

m 

(4.11) A f = ^A i V< ? ,(a:*); \ j9j (x*) = Vj = l,...,m. 

3 = 1 

(See e.g. [3].) We first prove by contradiction that if Slater's condition and the 
nondegeneracy assumption hold then Ao > 0. Suppose that Ao = and let J := 
{j G {1, . . . , to} : Xj > 0}; hence J is nonempty as A ^ 0. With x Q G K such 
that gj(xo) > (as Slater's condition holds, one such xo exists), let B(xo,p) ■= 
{z : \\z — xq\\ < p}. For p sufficiently small, B(xo,p) C K and gj(z) > for all 
z G B(xq, p) and every j = 1, . . . , m. Then by (|4.1ip and Ao = 0, 

m 

= J2 X 1 (^9j{x*),z-x*), Vz G B(xq, p), 

3 = 1 

which in turn implies (by nonnegativity of each term in the above sum) 

(Vgj(x*),z - x*) = 0, Vz G B(x Ql p) 1 j G J. 

But this clearly implies Vgj(x*) = for every j G J, in contradiction with the 
nondegeneracy assumption. Hence Ao > and by homogeneity, we may and will 
take A = 1. 

Therefore, letting Y := x* in (|4.6p . the polynomial X f(X) — f* can be 
written 

m 

ex-r = x>[<v<?i(<), 

3 = 1 



m 
3 = 1 



x*)g k (x)+ ^( x > 

x*)gk{x*) 

k=0 k=0,k^j 
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where we have used (14.61) with Y — x* and e = 0. Next, observe that : 



X»a jk (X,x*) G S 2 [X] [as G £ 2 [X, Y\] 
X h-> i> jk (X,x*)g k (x*) G S 2 [X] [as ^fc € S 2 [X, Y] and 9j (x* ) > 0] 
^j9j( x *) = i = l,...,m. 
And so, as A G Rip, 

(4.12) X^FX-f = ^o(X) + J2 A o( X )9 j (X), 

3 = 1 

for SOS polynomials (A^)^ C Y?[X] defined by 

m / m 

X^A (X) = ^A, £ il> jk (X,x')g k (x 

3=1 \k=0,k^j 
m 

X^AjiX) = ^Ajafc j = l,...,m. 

i=i 

Write every affine polynomial / G M.[X] as f T X + /q for some f G R™ and 
/o = /(0). If / is nonnegative on K then from (|4.12p , 

m 

/(x) = f t x - r + r +f = r + k + a po + J2 ^ w 9 3 po 

m 

i=i 

with A G £ 2 [X] (because f* + fo > 0) and so, the PP-BDR property holds for 
K with order d. By [TBI Theor. 2], K is SDr with the semidefmite representation 

621). □ 

We next show that the two sufficient conditions of strict convexity and SOS- 
convexity of Helton and Nie [6] in Theorem 14.21 both imply that Assumption 14.61 
holds and so Theorem 14.81 contains Theorem 14.21 as a special case. 

Corollary 4.9. Let K in be convex and both Assumvtion \2. 1\ and Slater's 

condition hold. Assume that either —gj is SOS-convex or —gj is convex on K and 
— V 2 <7j >~ on K n {x : gj(x) = 0}, for every j = 1, . . . ,m. Then Assumption^ 
holds and so Theorem \4-.8\ applies. 

Proof. By Lemma 12.51 for every j — 1, . . . , m, write 

(X,Y) i ► gj (X)-g(Y)-(V gj (Y),X-Y) 



(X-Y), 



[ [ V 2 g 3 {Y + s(X-Y)) dsdt) {X - Y) 
Jo Jo J i 



If —V 2 gj y on y G K with gj(y) — 0, then from the proof of [SJ Lemma 19], 
—Fj(x, y) y for all x, y G K with gj(y) = 0. In other words, —Fj(x, y) >; SI n on 
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Slj (defined in (|4.5|) ) for some 8 > 0. Therefore, by the matrix polynomial version 

of Putinar Positivstellensatz in 6, Theor. 29], 

(4.13) 

m m 

-F i {X,Y) = Y t a jk {X,Y)g k (X)+ £ $ jk {X, Y)g k {Y) + fa{X, Y) 9j (Y) 

k=0 fc=o,fc#j 

for some SOS matrix polynomials (dj k (X,Y)), ("ipj k (X,Y)) and some matrix poly- 
nomial ipj (X, Y) . 

On the other hand, if — gj is SOS-convex then by Lemma [2~4l —Fj(X, Y) is SOS 
and therefore (|4. 13|) also holds (take dj k = for all k ^ 0, \j}j k = for all k and 
Tpj = 0). But then 

gj (X)-g(Y)-(V 9j (Y),X-Y) = ((X — Y),Fj(X, Y)(X - Y)) 

m 

= -J2(( X - Y )^k(X,Y)(X-Y))g k (X) 

m 

- ((X -Y),$ jk (X,Y)(X -Y))g k (Y) 

k=0,k^ 3 

- {(X-Y) 1 ^(X,Y)(X-Y))g, J (Y) 

m 

k=0 
m 

^k (X, Y) g k (Y) - ^ (X, Y) g 3 (Y) 

k=0,k=£j 

for all X, Y and for some SOS polynomials djk, ipjk G M.[X, Y] and some polynomial 
ipj e R[X,F]. Equivalently, 

m 

(V 9j (Y),X-Y) = gj (X)- gj (Y)+J2^k(X,Y)g k (X) 

m 

+ Y ^k(X,Y)g k (Y)+^(X,Y) gj (Y) 

m m 

= 5>^pf,Y)< ?fc pO+ ^MX,Y)g k (Y) 

k=0 k=Q,k^j 

+ 1 p' j (X,Y)g j (Y) 

for some SOS polynomials cr'- k ,tpjk G S 2 [X, Y] and some polynomial ip'j € M.[X,Y]. 
In other words, Assumption 14.61 holds, which concludes the proof. □ 

Hence if each —gj is SOS-convex or convex on K with — V 2 <?j >- on K n {a; : 
gj{x) = 0}, one obtains a numerical scheme to obtain the parameter d in Thco- 
rem l4.8l as well as the semidefinite representation (|4.2[) of K. Solve the semidefinite 
programs (|4.8p with degree parameter dj. Eventually, pj — for every j = 1, . . . , m. 
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Example 3. Consider the convex set K in (|4.9p of Example [2] for which the defining 
polynomial g\ of K is not concave. We have seen that Assumption 14.61 holds (up to 
pi ~ 10~ n , close to machine precision) and max[di,d2] = 3. By Theorem 14.81 if pi would 
be exactly 0, the set 



(4.14) 



n := ^ 



M 3 (y) 
M 2 (gjy) 

Ly(Xi) 

yo 



h 
>0, 
= x l 
= 1 



1,2 
= 1,2 



would be a semidefinite representation of K. 

At least in practice, for every linear polynomial / G R[A], minimizing L y (f) over Q 
yields the desired optimal value /* := min^gK f(x), up to pi ~ — 10 _ . 

Indeed, let / G WL[X] be i T X for some vector f G K n . In minimizing / over K, one has 
f = Ai Vgi (x*) +A 2 Vg2{x*) for some A G R+, some x* G K with \igi(x*) = 0, i = 1, 2, and 
/* = Ai (Vgi(x*), x*) + \2(Vg2(x*),x*) = mm xe K.i T x. Let x be as in (I4.14p . arbitrary. 
Then 

2 

f T x-f =L y (f(X)-f) =J2*i L y({V9i(x*),X-x*)). 
If Ai > so that gi{x*) = 0, use (|4.12|) to obtain 

2 

M(V 5 i(x*),X - = L y (pi + A (A) + ^ A J (X)g 3 (X)) > pi, 

7=1 

because L y (Ao) > follows from A/a(y) >; 0, and L y (Ajgj) > 0, j — 1,2, follows from 
M 2 (sny), M 2 (p 2 y) b 0. If A 2 > so that 32(2:*) = 0, then from ftjfy 

L y {(yg 2 (x*),X- x*)) = L y (g 2 (X) - ((X - x"), V 2 g 2 (x*)(X - x*))) > 0, 

because L y {g2) > follows from M2(g2 y) b whereas the second term is nonnegative as 
{(X - x*), -V 2 g 2 {x*)(X - x*)) is SOS and M 3 (y) t 0. Hence f T i - /* > pi. On the 
other hand, from K C {x : (x,y) G Q}, one finally obtains the desired result 

f+pi < mm{f T x : (x,y) G fl} < f*. 



5. Conclusion 

As well-known, convexity is a highly desirable property in optimization. We have 
shown that it also has important specific consequences in polynomial optimization. 
For instance, for polynomial optimization problems with SOS-convex or strictly 
convex polynomial data, the basic SDP-relaxations of the moment approach [13] 
recognizes convexity and finite convergence occurs. Similarly, the set K has a 
semidefinite representation, explicit in terms of the defining polynomials (gj). 

The class of SOS-convex polynomials introduced in Helton and Nie [B] is partic- 
ularly interesting because the semidefinite constraint to handle in the semidefinite 
relaxation only involves the Hankel-like moment matrix which does not depend on 
the problem data! Hence one might envision a dedicated SDP solver that would 
take into account this peculiarity as Hankel-like or Toeplitz-like matrices enjoy very 
specific properties. Moreover, if restricted to this class of polynomials, Jensen's in- 
equality can be extended to linear functionals in the dual cone of SOS polynomials 
(hence not necessarily probability measures). 

Therefore, a topic of further research is to evaluate how large is the subclass of 
SOS-convex polynomials in the class of convex polynomials, and if possible, to also 
provide simple sufficient conditions for SOS-convexity. 
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